arXiv:1502.01881vl [cond-mat.quant-gas] 6 Feb 2015 


Scattering of matter-waves in spatially inhomogeneous environments 


F. Tsitoura, 1 P. Kruger, 2 P. G. Kevrekidis, 3,4 and D. J. Frantzeskakis 1 

/ Department of Physics, University of Athens, Panepistimiopolis, Zografos, Athens 15784, Greece 
2 Midlands Ultracold Atom Research Centre, School of Physics & Astronomy, The University of Nottingham, Nottingham, UK 
3 Department of Mathematics and Statistics, University of Massachusetts, Amherst MA 01003-4515, USA 
4 Center for Nonlinear Studies and Theoretical Division, 

Los Alamos National Laboratory, Los Alamos, NM 87544, USA 

We study scattering of quasi one-dimensional matter-waves at an interface of two spatial domains, one with re¬ 
pulsive and one with attractive interatomic interactions. It is shown that the incidence of a Gaussian wavepacket 
from the repulsive to the attractive region gives rise to generation of a soliton train. More specifically, the num¬ 
ber of emergent solitons can be controlled e.g. by the variation of the amplitude or the width of the incoming 
wavepacket. Furthermore, we study the reflectivity of a soliton incident from the attractive region to the repul¬ 
sive one. We find the reflection coefficient numerically and employ analytical methods, that treat the soliton as 
a particle (for moderate and large amplitudes) or a quasi-linear wavepacket (for small amplitudes), to determine 
the critical soliton momentum - as function of the soliton amplitude - for which total reflection is observed. 
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I. INTRODUCTION 

For almost two decades, the study of nonlinear phenom¬ 
ena occurring in atomic Bose-Einstein condensates has expe¬ 
rienced an enormous increase of interest mu. A prominent 
example, in the quasi one-dimensional (ID) setting, is the ex¬ 
perimental observation of robust matter-wave solitons of the 
bright H and dark 01 type, and the study of their properties 
(see, e.g., the reviews 0 and 0 for bright and dark solitons, 
respectively). Such coherent nonlinear excitations of BECs 
are also interesting from the viewpoint of potential applica¬ 
tions, ranging from coherent matter-wave optics to precision 
measurements and quantum information processing. Indeed, 
the formal similarities between nonlinear and matter-wave op¬ 
tics 0 indicate that coherent matter-waves may in princi¬ 
ple be controlled similarly to their optical siblings in optical 
fibers, waveguides, photonic crystals, and so on 0. 

In that respect, it is not surprising that there exist many 
works devoted to the manipulation of matter waves. Among 
various techniques that have been proposed, an experimen¬ 
tally tractable one refers to engineering the “environment” of 
the matter-wave, by magnetically 0 or optically El induced 
Feshbach resonances, which makes it possible to control the 
effective nonlinearity in the condensate. The application of 
such a “Feshbach resonance management” (FRM) technique 
O in the temporal domain was used for the realization of 
matter-wave bright solitons by switching the interatomic in¬ 
teractions from repulsive to attractive 0; it was also pro¬ 
posed as a means to stabilize attractive higher-dimensional 
BECs against collapse [12, O] and to create robust quasi-ID 
matter-wave breathers fill fl4ll . On the other hand, the FRM 
technique in the spatial domain, which gives rise to the so- 
called “collisionally inhomogeneous condensates’ ’ E! with 
a spatially modulated nonlinearity, has also been extensively 
studied. In particular, novel phenomena and a variety of ap¬ 
plications have been proposed in this context, including the 
adiabatic compression of matter-waves El IH, Bloch oscil¬ 
lations of matter-wave solitons El, atomic soliton emission 
and atom lasers El, enhancement of transmittivity of matter 


waves through barriers [18], formation of stable condensates 
exhibiting both attractive and repulsive interatomic interac¬ 
tions [ |19fl , solitons in combined linear and nolinear potentials 
[20||, generation of_solitons [21] and vortex rings 1221 . control 
of Faraday waves l23ll . and many others. A detailed recent re¬ 
view of such inhomogeneously nonlinear settings, especially 
in the context of periodic (i.e., nonlinear lattice) variations can 
be found in JH. 

In this work, we study the scattering of matter-waves in a 
collisionally inhomogeneous environment. In particular, we 
consider a quasi-ID setting (whereby matter waves are ori¬ 
ented along the x-direction) and assume that the scattering 
length a is piecewise constant for and i>0, taking 

respectively the values — a\ < 0 and 02 > 0, and changes 
sign at x = 0. In other words, we assume that the normalized 
scattering length a(x) takes the form: 





, (1) 


where W is the spatial scale over which the transition from 
the asymptotic values a\ and a 2 takes place. For the above 
setting, and in the framework of the mean-field approxima¬ 
tion, we will investigate two different scattering processes; a 
description of our considerations and the organization of the 
paper are as follows. 

First, in Section II, we study the incidence of a nearly linear 
(Gaussian) wavepacket from the repulsive region (x > 0) to 
the attractive region (x < 0), and demonstrate the generation 
of a train of bright solitons. By numerically integrating the 
pertinent Gross-Pitaevskii (GP) equation, we determine the 
number of created solitons as functions of the initial data (am¬ 
plitude, width and momentum of the incident wavepacket), as 
well as the difference of the values of the scattering length. 

In Section III, we study the reflectivity of a bright soliton 
from the scattering length interface; the soliton is assumed to 
exist and travel from the attractive region (x < 0) towards the 
repulsive region (x > 0). We find numerically the reflection 
coefficient as a function of the soliton momentum and ampli¬ 
tude, and find that it has a step-like dependence on momentum 
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for sufficiently weak solitons. In the case of total reflection, 
we use an analytical approximation (treating the soliton as a 
particle) and find the equation of motion for the soliton cen¬ 
ter. This equation is used to determine the critical value of 
momentum below which total reflection occurs, which turns 
out to depend linearly on the soliton amplitude. Addition¬ 
ally, for extremely weak solitons, employing results from lin¬ 
ear quantum mechanics iBl , we also find a (different) linear 
dependence of the critical momentum on the soliton ampli¬ 
tude. Both analytical estimates, for weak and strong solitons 
are found to be in very good agreement with the numerical 
results, with the latter also encompassing a transition region 
between the two regimes. 

Finally, Section IV concludes our findings and presents a 
number of directions for future study. 


II. REFLECTIVITY OF A GAUSSIAN WAVEPACKET 
FROM THE SCATTERING LENGTH INTERFACE 

A. Model and creation of a soliton train 


Our considerations start from the following Gross- 
Pitaevskii (GP) equation, which describes a quasi-ID BEC 
oriented along the x-axis 0111: 



h 2 d 

2m dx 2 


2Huj±a\'i’\ 2 'i!. 


( 2 ) 


Normalizing time and space in Eq. 0. as t —> u±t and 
x —> x/a± (where a± = ( H/mcj± j 1 / 2 is the transverse har¬ 
monic oscillator length), as well as the density as \i/j\ 2 —> 
2cti|u| 2 , we cast Eq. (0 into the following dimensionless 
form: 


.du 1 d 2 u . .. ,2 

'm + a(l)w " = °’ <4) 

where the function a(x) is given by Eq. 0. For our analytical 
and numerical considerations below, we will use the values 
< 22/01 = 0.95, and W = 0.01 corresponding to an abrupt, 
step-like transition. 

It is relevant to point out here that, generally, Eq. 0 as 
well as its variants attempting to more adequately capture 
transverse degrees of freedom (see the quasi-ID models of 
Refs. 0 and 130]), suggest that the variations/modulations 
in transverse trapping strength can be used in a way equivalent 
to longitudinal variations of the scattering length. This idea 
has been used even in a quantitative fashion, e.g., to explain 
the phenomenology of the formation of Faraday wave patterns 
- cf. Ref. ll3lll (and relevant work in Ref. El ). Nevertheless, 
this type of consideration is not applicable in the present set¬ 
ting, given the sign changing nature of the nonlinearity. 

We now assume that a Gaussian wavepacket of amplitude 
Uq and width l, initially located at x = xq > 0 (i.e., in the 
repulsive region), moves towards the attractive region. The 
specific form of the wavepacket, which is used as an initial 
condition for Eq. 0 in our simulations, is: 


Here T<(x,f) is the mean-field order parameter, m is the 
atomic mass, w_l is the transverse confining frequency, and 
a the s-wave scattering length [a > 0 (a < 0) corresponds to 
repulsive (attractive) interatomic interactions]. Considering a 
situation where proper spatially dependent fields close to Fes- 
hbach resonances are employed, we assume that the scattering 
length a is piecewise constant for x < 0 and x > 0, taking the 
form 

a{x) = (l/2)[(a 2 — ai) + (ai + a 2 ) tanh(x/kF)], (3) 

where W is the spatial scale over which the transition from 
the asymptotic value — ai < 0 (for x/W —> — 00 ) to a 2 > 0 
(for x/W —»■ + 00 ) takes place. 

A strategy for developing a corresponding experimental im¬ 
plementation can be based on the interaction tunability of spe¬ 
cific atomic species by applying external magnetic fields. For 
example, for cesium the s-wave scattering length a changes 
sign through a zero-crossing at an external field strength of 
17 G [26]]. Confining cesium atoms in an elongated trapping 
potential near the surface of an atom chip [27] will allow for 
appropriate local engineering of a to form steps of varying 
widths W, where the atom-surface separation sets a scale for 
achievable minimum step widths. The trapping potential can 
be formed optically, possibly also by a suitable combination 
of optical and magnetic fields, whereby care has to be taken 
as the magnetic field will influence both the external potential 
and the scattering length profile a(x)\ see e.g. the relevant 
discussion of 128], 


u(x, 0) = Uo exp j exp (-iKx), (5) 

where K is the initial momentum of the wavepacket. No¬ 
tice that this form of the wavepacket approximates the ground 
state profile in the case of relatively small atom numbers (cor¬ 
responding to a weak nonlinearity). 

Using the parameter values xo = 20, l = 10, and K = 1.5 
(as well as Uo = 1, W = 0.01 and a 2 /ai = 0.95), we depict 
the corresponding configuration in the top panel of Fig. □ In 
the middle and bottom panels of the same figure, we show the 
subsequent dynamics: it is observed that the wavepacket is 
transmitted through the discontinuity of the scattering length 
at x w 0 and, after entering the region with attractive interac¬ 
tions, it transforms into a train of bright solitons. Notice that 
the soliton generation process is such that each generated soli¬ 
ton is larger than the one that will be generated at a later time. 
This is due to the fact that once a portion of the condensate 
enters the attractive side and is self-organized into a soliton, 
the number of atoms of the wavepacket on the repulsive side is 
decreased and, thus, a smaller soliton will be generated next. 
It is interesting to observe that the ratio of the velocities of two 
adjacent solitons in the train is constant; as a result, for each 
certain time instance, the distance between adjacent solitons is 
the same. Note that the results on the generation - and charac¬ 
teristics - of the atomic soliton train described above are remi¬ 
niscent to the ones found in Ref. El , but by means of a some¬ 
what different physical mechanism: in that work, the soliton 
train was produced via a sufficiently deep spatially-dependent 
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FIG. 1: (Color online) Motion of a Gaussian wavepacket from the 
region with repulsive interactions to the region with attractive inter¬ 
actions. Top and middle panels show, respectively, the density pro¬ 
files of the wavefunction at at t = 0 and t = 200 [solid (blue) lines]; 
the normalized scattering length profile for W = 0.01 is also de¬ 
picted [dashed (red) line]. The bottom panel represents a contour 
plot showing the evolution of the density; middle and bottom panels 
clearly show the creation of a bright soliton train. Parameter values 
are: Uo = 1, xo = 20, l = 10, K = 1.5, and 02/01 = 0.95. 


nonlinearity which acted on a trapped Gaussian wavepacket 
(existing between a region of vanishing and that of negative 
scattering length). The depth of the (abrupt) negative step was 
found to control the number of emitted solitary waves. 

We find that the number of the created solitons, N s , de¬ 
pends on the momentum K, the amplitude Uq and the width 
l of the Gaussian wavepacket, as well as the height of the in¬ 
terface 02 / 01 . Results pertaining to the count of the soliton 
number are shown in Fig. [2] larger initial amplitudes and/or 
widths of the wavepacket result in a larger number of solitons. 
On the other hand, increasing the initial momentum k and/or 



FIG. 2: (Color online) The top and bottom panels show, respectively, 
the number of solitons, N s , observed at t = 200, as a function of the 
initial amplitude Uo (for fixed l = 10) or as a function of the initial 
wavepacket width l (for fixed Uo = 1) of the Gaussian wavepacket. 
The other parameter values are W = 0.01, ai = 1, <22 = 0.95, 
K = 1.5, and xo = 20. 


the height 02/01 of the interface the number of solitons seems 
to have a weaker effect on the process; for the particular ex¬ 
ample shown in Fig. [Q the number of solitons is 7 at time 
t = 200. Here, we should note that for the counting of the 
number of solitons, we have included only solitons of ampli¬ 
tudes at least 10 % of the first created soliton. 


III. REFLECTIVITY OF A SOLITON FROM THE 
SCATTERING LENGTH INTERFACE 


Next, we consider the reflectivity of a bright soliton at the 
scattering length interface. We assume, in particular, that a 
bright soliton moves from the attractive (x < 0 ) to the repul¬ 
sive region (x > 0 ) and is, thus, scattered at the interface, at 
x = 0 , caused by the change of the sign of the nonlinearity. 
This dynamical scenario is, effectively, complementary to the 
one studied in Section II. 

The bright soliton propagating in the attractive region has 
the form; 

u(x,t) = rjsech [77 (x — xo(f))] exp (i (kx — uit )), ( 6 ) 

where rj, k, Xq and 00 respectively denote the amplitude, ve¬ 
locity, initial position and frequency of the soliton. Then, we 
numerically integrate Eq. © with the initial condition taken 
as 


u(x, 0) = Tysech [rj (x — xo(0))] exp ( ikx ) exp (/</>), (7) 

and observe the dynamics of the scattering process. Typical 
outcomes are shown in Fig. [3j in all cases, we fix the initial 
soliton momentum, at k = 0.7, and vary the amplitude 77 . 
We observe that if the soliton amplitude is sufficiently small 
(large) then total transmission (reflection) is found -cf. top 
(bottom) panel of the figure for ij = 0.2 (rj = 1). On the other 
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FIG. 4: (Color online) The reflection coefficient R as a function of 
the initial soliton momentum k, for various values of large soliton 
amplitudes tj. The inset shows cases corresponding to weak solitons. 


FIG. 3: (Color online) Contour plots showing the evolution of the 
density of a bright soliton of initial velocity k = 0.7, scattered at the 
interface between the attractive and repulsive region, at x = 0 [de¬ 
picted by the dashed (white) line]. From top to bottom, the amplitude 
of the soliton is 77 = 0.2, rj = 0.6 and r/ = 1. Top (bottom) panel 
corresponds to total transmission (reflection); middle panel shows 
partial reflection. 


hand, for a moderate value of 77 (e.g., rj = 0.6 -cf. middle 
panel) the soliton is partially transmitted and reflected. 

The soliton reflectivity can be calculated numerically upon 
determining the reflection coefficient R, defined as the num¬ 
ber of atoms remaining in the x < 0 (attractive) region over 
the number of atoms of the incident soliton. Taking into re¬ 
gard that the latter is given by J | u(x, 0 ) \ 2 dx = 277 , we can 
express R as: 


1 f u 

R = — \u(x,U)\ 2 dx. ( 8 ) 

2?7 J -00 

Here, f* is a time sufficiently large such that the reflected and 
transmitted parts of the soliton are spatially well separated; 
this separation is set by a spatial region of extent Ax « /cf* 
around x = 0 and, accordingly, f, is appropriately chosen for 
each individual numerical experiment. 

Figure [4] shows the reflection coefficient as a function of 
the initial soliton momentum (for 0 < k < 10 ) and vari¬ 
ous values of the soliton amplitude 77 . We observe that when 
77 is increased, the respective reflection coefficient curves 
drift towards larger momentum values and the curves become 
smoother: the transition from total reflection to total trans¬ 
mission becomes less sharp. This means that the interval of 
momenta for which partial transmission and reflection occur 
(as in the middle panel of Fig. 0 increases with increasing 
soliton amplitude. 

From the above discussion, it is obvious that the soliton 
keeps its particle-like character only in the case where it is 
totally reflected (cf. bottom panel of Fig. 0: in the cases of 
total or partial transmission, the soliton is not supported in the 
repulsive regime and it is eventually destroyed. We can thus 
adapt the particle picture for the soliton dynamics in the total 


reflection regime, and describe analytically the soliton trajec¬ 
tory and its reflectivity properties. Our approach based on 
the center of mass (defined below) extends the corresponding 
considerations of Ref. |33l l, where a similar methodology was 
developed for the case of a linear step potential. 

We start with the soliton’s center of mass, given by: 

/ OO 

x\u(x,t)\ 2 dx, (9) 

-OO 


which is connected with the soliton momentum P = 
f^° oo (uu*—u*u x )dx through the equation dx/dt = P. 
Then, differentiating the latter expression with respect to f, 
and using Eq. <0, it is straightforward to derive the following 
equation of motion for x: 


d 2 x 

~dP 


1 

2 



da{x) 

dx 


| u(x, t)\ A dx. 


( 10 ) 


The integral on the right-hand side of Eq. (flOb can be cal¬ 
culated in an analytical form, upon approximating a(x) [cf. 
Eq. 0] by a Heavyside function in the limiting case where 
W -C I/77. Then, da(x)/dx is approximated by a delta func¬ 
tion, and integrating the right-hand side of Eq. (flOb . we end 
up with the following result: 


“ 4( 1 + S) ,W<,x ° <i)) ' (H) 


Then, taking into regard that the soliton center is connected 
with the center of mass through the equation x - 277 x 0 , we 
can express Eq. CEB as follows: 

d 2 x 0 dV efi 

dt 2 dx o’ ( ’ 

where the effective potential V e g is given by: 


FetrE-t-o) — 


x 


1 

12 





2 


[3tanh ( 77 x 0 ) — tanh 3 ( 77 x 0 )] ■ 


(13) 
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FIG. 5: (Color online) A sketch of the effective potential 14g (blue 
line) as a function of xo■ Shown also is the soliton, initially located 
at *o(0) [left (green) curve], far from the scattering length interface, 
and in close proximity to the interface [right (red) curve], where the 
location of its center is xo(| Ax|). 


Equation (fl2l > shows that the soliton can be regarded as a New¬ 
tonian unit-mass particle, which evolves in the presence of the 
effective potential V e g', the latter, has a shape of a step-like 
barrier, as depicted in Fig. [5] Thus, according to this par¬ 
ticle picture, the soliton will be totally reflected if its initial 
energy E s is less than the “height” of the barrier. Since the 
soliton is expected to interact with the effective potential only 
through its exponential leading tail, the soliton center is an¬ 
ticipated to never reach the interface at x = 0 , but it will 
approach it only up to a distance roughly equal to the half¬ 
width at half-maximum (HWHM) of the soliton; the above 
situation is schematically illustrated in Fig. [5] Thus, tak¬ 
ing into regard that the soliton’s HWHM, denoted by Ax, 
is connected with the inverse width rj through the equation 
Ax = ln(l + y/2 )/r), we can find that the relevant barrier 
height is given by I4ff(Ax) — I4ff(xo(0)), where xo(0) is the 
initial soliton position. According to the above arguments, the 
soliton will be totally reflected if the initial soliton energy is 
less than (or equal to) the effective barrier height, namely: 

E s = h 2 + V eS (x 0 {0)) < T4ff(Ax) (14) 

We have numerically checked the validity of this analysis 
by comparing, at first, the numerically obtained soliton trajec¬ 
tory [by means of direct numerical integration of Eq. Q in the 
case of total reflection] with the approximate analytical result 
ofEq. O. A typical example, corresponding to a soliton am¬ 
plitude 77 = 0.4 and momentum k = 0.1, is shown in Fig. [ 6 ] 
There, the numerical result is displayed in the form of a con¬ 
tour plot for the evolution of the soliton density, as well as the 
analytical result of Eq. (IT2l > - cf. dashed line in the figure. 
Note that similar results were obtained for soliton amplitudes 
0.2 < 77 < 2. It can be seen that the dashed line follows with 
a fairly good accuracy the evolution of the soliton center. The 
slight discrepancy observed can be explained as follows: the 
tail of the bright soliton, in case of total reflection (cf. Figs.0 
and [ 6 j, interacts with the interface and enters the repulsive 
area, and eventually comes back to the attractive region. This 
effect, which cannot be explained via the particle approach, 
causes a slight shift in the soliton trajectory. Thus, the trajec- 
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FIG. 6 : (Color online) Contour plot showing the space-time evolu¬ 
tion of the motion of the soliton center. The bright soliton is initially 
placed at xo(0) = —10, while its amplitude and momentum are 
77 = 0.4 and k = 0.1. The dashed (black) curve represents the an¬ 
alytical result of Eq. while the horizontal (white) line depicts 
x — 0 . 


tory obtained from Eq. (ITOl) has naturally a slight discrepancy 
for any soliton amplitude 77 . 

Next, employing Eq. (fl4l) . it is possible to derive analyti¬ 
cally the critical value of the initial momentum fc cr [when the 
equality in Eq. ( fl4l > holds], for which total reflection occurs, 
as a function of the soliton amplitude and the parameters char¬ 
acterizing the scattering length profile. The result is: 


k 


cr 



1/2 


(i + c) 


V, 


(15) 


where constant C = (l/2)[ 3 tanh^Ax) — tanh 3 ( 7 ;Ax)] and 
77 AX = — ln(l + -\/2) ~ —0.88. Note that Eq. (IT5l) suggests a 
linear dependence of k CI on 77 , which is confirmed by our nu¬ 
merical simulations. Indeed, as shown in Fig. [7J for soli tons 
of sufficiently large amplitudes, i.e., for 77 > 0 . 2 , this ana¬ 
lytical prediction [depicted by the solid (green) straight line] 
is in an excellent agreement with the numerical result for fc cr 
[depicted by the (red) dots]. Notice that the numerically ob¬ 
tained values for k CI are calculated so that the respective re¬ 
flection coefficient values become less than unity by a factor 
of 10 ~ 3 ; however, we note here that the results presented are 
only weakly sensitive to the selection of the particular thresh¬ 
old. 

For weaker solitons it is expected that our analytical ap¬ 
proximations described above should be less accurate: this is 
due to the fact that for small values of 77 , the nonlinearity be¬ 
comes extremely weak, and thus a linear description of the 
problem would be more appropriate. In such a case, the soli¬ 
ton can be treated as a linear wavepacket, which is scattered 
from an effective step barrier; the latter, is basically formed 
by the step-like change of the scattering length profile. Then, 
the reflection coefficient can be approximated from the corre¬ 
sponding linear problem I 25ll as follows: 


4 y/(E - V 0 )E 

{VE + Ve =%) 2 ’ 


( 16 ) 
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FIG. 7: (Color online) The critical momentum k CI , below which total 
reflection occurs, as function of the soliton amplitude r/. The (red) 
dots depict the numerical result, while the left (blue) and right (green) 
straight lines correspond to the analytical predictions of Eqs. d!9t and 
GS>, respectively. The parameter values are a\ = 1, a ,2 = 0.95 and 
W = 0.01. 

where E and Vo denote, respectively, the energy of the 
wavepacket and the height of the effective potential barrier. 
Notice that Eq. ( fT 6 l > stands for plane waves; however, it can 
still provide a reasonable approximation as long as the soli¬ 
ton width ? 7 _1 is sufficiently large, i.e., for sufficiently weak 
solitons. In our case, the soliton energy is given by (see, e.g.. 
Ref. i]): 

E = rjk 2 - i? 7 3 , (17) 

while the strength of the effective barrier potential is given by: 

Vo = \j a(x)\u\ A dx = i ^ - lj rf. (18) 

Then, total reflection, i.e., R = 1 in Eq. (IT 6 l) . occurs for E = 
Vb; the latter equation leads to the following result for the 
critical momentum 

= (£•) v - (19) 

The above approximate analytical result, which is relevant to 
weak solitons, also shows a linear dependence of k CI on 77 
and is in a very good agreement with the numerical results, as 
shown in Fig.|7]ior r) < 0.1. 

In summary, we capture the regime of small 77 by means of 
the linear/wave theory, and the regime of large 77 by our soliton 
particle theory, while between the two we interpolate via the 
use of numerical computations as shown in Fig. [7] 

IV. CONCLUSIONS 

In this work, we studied the scattering of quasi-ID matter- 
waves in a spatially inhomogeneous environment, character¬ 
ized by a piecewise constant profile of the scattering length a, 
such that a = —a\ < 0 for x < 0 , a = a 2 > 0 for x > 0 , 
and a changes sign at x = 0. This way, in the region x < 0 
(x > 0) the interatomic interactions are attractive (repulsive). 
We investigated two different dynamical scenarios: 


(i) the scattering of a quasi-linear (Gaussian) wavepacket 
at the scattering length interface, with the wavepacket 
traveling from the repulsive to the attractive region, and 

(ii) the scattering of a matter-wave bright soliton at the scat¬ 
tering length interface, with the soliton traveling from 
the attractive to the repulsive region. 

In case (i), we found that when the wavepacket enters 
the attractive region it evolves into a train of bright soli¬ 
tons. The soliton train is such that each generated soliton 
is larger than the one that will be generated at later times, 
while the distance between adjacent solitons is the same. We 
counted the number of the created solitons, as a function of 
the wavepacket’s initial characteristics (momentum, ampli¬ 
tude, and width) and the height of the nonlinearity interface 
02 / 01 , and found that larger initial amplitudes and/or widths 
of the initial wavepacket result in a larger number of solitons. 

For case (ii), we found that the incidence of the soliton at 
the scattering length interface generally leads to total trans¬ 
mission, total reflection, or partial transmission/reflection. 
The reflection coefficient was determined numerically as a 
function of the initial soliton momentum, for different soli¬ 
ton amplitudes. For sufficiently weak solitons, we found an 
almost abrupt change from total transmission to total reflec¬ 
tion, effectively associated with the linear phenomenology in 
a step potential. For stronger solitons, the reflection coeffi¬ 
cient featured a smoother dependence on the momentum. 

We also developed analytical approximations - that treated 
the solitons as particles (for large amplitudes) or linear 
wavepackets (for small amplitudes) - to determine the crit¬ 
ical value of soliton momentum, k cr , below which total re¬ 
flection occurs. We found that k CT depends linearly on the 
soliton amplitude, but with different slopes in the purely non¬ 
linear and the quasi-linear regimes. Numerically, we find a 
smooth crossover between these two regimes, which can be 
interpreted as a gradual continuous change of the soliton from 
being dominated by wave-like to particle-like properties. Our 
analytical predictions were found to be in very good agree¬ 
ment with the corresponding numerical results. 

There are numerous directions that may be worth consid¬ 
ering for future efforts. One of these is to consider the pos¬ 
sibility of multiple steps and their interplay. Another is to 
examine the interplay of the nonlinear step with an external 
linear potential or with a non-trivial background (e.g. on the 
repulsively interacting side which can support such a back¬ 
ground) and to explore the dynamics of incident wavepackets 
in such settings. Potentially, probing the soliton dynamics in 
such configurations could be utilized towards retrieving quan¬ 
titative information about the nature of linear and/or nonlinear 
unknown potentials. 

From a more rigorous mathematical perspective, it will be 
interesting to attempt to connect the present setting to the 
extensive developments in treating integrable problems with 
suitable boundary conditions (e.g. on the half line), as de¬ 
tailed e.g., in Ref. ll34fl . A way to make this connection may 
be to consider the GP equation e.g. solely on the attractive 
domain with a boundary condition inferred by the incidence 
of the Gaussian wavepacket at x = 0 (i.e., a Gaussian in time 
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boundary condition). A potential by-product of such a for¬ 
mulation might be the identification of the number of solitary 
waves that will emerge, as a function of the properties of this 
effective (and localized in time) boundary drive. 

Finally, it would be of particular interest to extend con¬ 
siderations to the two- or higher-dimensional setting. There, 
understanding the properties of the formed solitons, e.g., on 
the attractive interaction “domain”, taking into consideration 
the collapse feature that arises in the critical or super-critical 
higher-dimensional case |3j|, would be especially relevant. 
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